Distribution of conductance for Anderson insulators: A theory with a single parameter 
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We obtain an analytic expression for the full distribution of conductance for a strongly disordered 
three dimensional conductor within a perturbative approach based on transfer matrix formulation. 
Our results confirm the numerical evidence that the log-normal limit of the distribution is not reached 
even in the deeply insulating regime. We show that the variance of the logarithm of the conductance 
scales as a fractional power of the mean, while the skewness changes sign as one approaches the 
Anderson metal-insulator transition from the deeply insulating limit, all described as a function of 
a single parameter. The approach suggests a possible single parameter description of the Anderson 
transition that takes into account the full non-trivial distribution of conductance. 

PACS numbers: 73.23.-b, 71.30., 72.10. -d 



Quantum fluctuations have been intensely studied in 
recent years, but in many cases a fundamental under- 
standing of their effects on physical observables remain 
poorly understood. In particular the effects of large 
mesoscopic fluctuations on quantum phase transitions, 
both in interacting electron systems^ and in disordered 
non-interacting models^, have not been studied system- 
atically. For a non-interacting system, while the distri- 
bution P(g) of the dimensionless conductance g in the 
metallic weak disorder limit is well understood^, even 
a qualitative understanding of P(g) is lacking in the 
strongly disordered insulating regime in three dimen- 
sions (3D) where numerical data show large deviations 4 
from an expected log-normal distribution. The more fun- 
damental question of how this non-trivial distribution 
changes as one decreases the disorder from the deeply 
insulating regime toward the Anderson metal-insulator 
transition point has remained largely unexplored. This 
is primarily due to the lack of appropriate theoretical 
tools to consider such distribution functions analytically. 
Conventional field theory framework* which relies on a 
small e expansion in 2 + e dimensions gives results that 
do not agree even qualitatively, in the e — > 1 limit, with 
numerical results in 3D£. 

In this work, we obtain the full conductance distri- 
bution in 3D in the strong disorder regime analytically 
within a perturbative approach, and show that the results 
are consistent with available numerical data. Our model 
involves a single disorder parameter T = £ /4L Z where £ is 
the localization length and L z is the length of the conduc- 
tor. A second parameter 7 = £/8L, where L is the cross- 
sectional dimension, becomes independent only if arbi- 
trary geometrical shapes are considered. We show first 
of all that the interaction between different channels of 
a conductor in 3D remains important even in the deeply 
insulating regime; the result is that the 3D distribution 
is never log-normal. Instead, we find that the variance 
of the logarithm of conductance scales approximately as 
2/3rd power of the mean 7 -. In addition, we also explicitly 
evaluate the third cumulant K3 = (Qxig— (lag)) 3 ) which 
describes the asymmetry of the distribution. We find 
that in the deeply insulating regime K3 is positive (longer 
tails toward larger conductances). It then decreases with 



decreasing disorder according to K3 ~ (—lug), going 
through zero well before the Anderson metal-insulator 
transition points. As one decreases the disorder further 
toward the critical point, K3 becomes negative, describ- 
ing increasing asymmetry in the opposite direction^. We 
emphasize that all of this is described as a function of a 
single parameter that fixes the mean value. Our method 
therefore allows us for the first time to explore in quan- 
titative detail how the conductance distribution changes 
as one approaches the Anderson transition point starting 
from the deeply insulating regime, and how it can still 
be described analytically within a one parameter theory. 

Our formulation is based on the transfer matrix 
framework developed originally for transport in quasi 
one-dimension (Q1D) where the transverse length of 
a conducting wire is less than the localization length. 
In Q1D, the Dorokhov-Mello-Pereyra-Kumar (DMPK) 
equation 1 ^ has been enormously successful in describ- 
ing the details of the distribution of the transmission 
eigenvalues 11 . Exploiting the Landauer formula 1 ^ to con- 
nect the distribution of g with the distribution of the 
transmission eigenvalues, the DMPK equation has been 
used to obtain a variety of novel features in the distri- 
bution of conductances in QID 1 ^. A generalization of 
the Q1D DMPK equation, claimed to be valid in 3D, has 
been proposed in [14]]. By solving this so called Gener- 
alized DMPK (GDMPK) equation numerically and com- 
paring the results with those from direct numerical solu- 
tion of the tight binding Anderson model, it has recently 
been shown-^ that the GDMPK not only incorporates 
the effects of dimensionality correctly, but that it also 
describes the full distribution of the transmission levels 
quantitatively in the insulating as well as near the critical 
regime in 3D. However, analytic solutions of the equation 
in the strongly disordered regime have been obtained only 
within a very approximate saddle point scheme, where 
the interaction between the transmission eigenvalues are 
entirely ignored 1 ^. These approximate solutions do show 
deviations from the Q1D behavior in the right direction, 
but they fail to describe correctly how e.g. the variance 
or the skewness of the distribution changes with disorder. 

Here we solve the GDMPK equation analytically, in- 
cluding for the first time the interaction between the 
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eigenvalues, within a novel perturbative approach. We 
find that the density of the transmission eigenvalues in 
the insulating regime is a constant, with an exponential 
gap at the origin that increases with increasing disor- 
der. This is fundamentally different from the approxi- 
mate solutions where interaction between eigenvalues are 
neglected. The resulting analytic expression for the full 
P(ln g) as a function of disorder now agrees quantita- 
tively with available numerical data. The results suggest 
a possible single parameter description of the Anderson 
transition starting from the insulating side and taking 
into account the full distribution of conductances even 
when the average (or the most probable value) of the 
conductance is no longer a meaningful representation of 
a highly non-trivial asymmetric distribution. 

For an iV-channel disordered conductor of fixed cross 
section L 2 , the GDMPK equation 14 describes the evolu- 
tion with length L z of the joint probability distribution 
PL z (yi) of the N transmission eigenvalues xf. 
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with the initial condition p(x, t = 0) = 5(x), where 
t = L z /l with I being the mean free path, f{xi,Xj) = 
sinh 2 Xj — sinh 2 Xi and we have kept A = 1 as a free 
parameter for later convenience. Here 7^ = 2Kij/Ku, 
where K^ is a phenomenological matrix defined in terms 
of certain eigenvector correlations that can be explicitly 
evaluated numerically. The Q1D DMPK equation is re- 
covered when 7ii_= 1 (we only consider orthogonal sym- 
metry). In Ref 15] it was shown that only two parame- 
ters, Kn and K12, are enough to model the entire matrix 
Kij in the insulating as well as the critical regimes, as 
proposed in Ref Disorder is then characterized by 
the parameter Y = ljK\\L z . In the insulating side, one 
can interpret £ = Al/Ku as the localization length, so 
T = £/4L z . The other parameter 712 = 7 = £/8£, so 
that Y/j = 2L/L Z depends on geometry. The entire dis- 
tribution, for a cubic system, is therefore characterized 
by a single disorder parameter Y (with 7 = Y/2). 

We note that it is through the L dependence of the 
eigenvector correlations Ku and K\2 that the GDMPK 
'knows' about the dimensionality of the system. In 3D, 
Ku ~ 1/LP (1/L 2 ) in the insulating (metallic) regime; 
the quantity K\\ = lim^^oo Ku(L) is zero in the metal- 
lic regime as well as at the critical point, but is finite 
for insulators 1 ^-. In 2D on the other hand, Ku is always 
finite. Thus the particular L dependence of K n, and 
therefore of our parameter 7, not only reflects the proper 
dimensionality, but also contains information about the 
critical point. 

Once an analytic solution of the GDMPK equation for 
p(x) is available, the full distribution of conductances 



P(g) can be obtained via the Landauer formula 12,17 

P(9) « J J] dx aP (x)5 - sech 2 x^ . (2) 

We can therefore expect that Eq. |T]) can be used as 
the starting point for studying the distribution of con- 
ductances at large disorder. Of course the DMPK itself 
breaks down for £ < I; however this happens far from the 
transition region and we do not consider such extreme 
disorder. 

We first briefly outline our method used to solve Eq. (TT]) 
in the insulating regime. Following Ref. [H[ we use a 
factorization in terms of f = e _S1 / 2 : 



p(x,t)= ; limC(x)G A r(x;t|y)r 1 (y) 
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nl sinh2iEi| 2 
Giv(x,t]y). (3) 
i I sum 2y, 1 2 



The GDMPK transforms into an equation for the N- 
particle Greens function Gn- —dG^jdt = HGn, with 
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where 7 = £/4L The initial condition is G(x; t = 0|y) = 
W.^Tr(y)^( x ~~ y) where ir(y) refer to a symmetrized 
permutation. This maps the problem onto a set of N 
interacting bosons evolving in imaginary time in ID with 
delta function initial conditions 19 . 

In the strongly disordered limit the interaction 
strength 7 <C 1, and we exploit this small parame- 
ter to develop a perturbation theory to evaluate the 
non-equilibrium Greens function using standard Keldysh 
techniques 2 ^. However, the small y singular behavior 
in the denominator of Eq. (|3|) demands that Gat(x, t\y) 
must go as ~ £(y) in the small y limit in order to recover 
a well-defined p. To extract this non-analytic behavior of 
Gn from the diagrammatic expansion, we must sum the 
expansion in a particular way. First we treat A = (A — 2) 
in Eq. (|4]) as a small parameter in which we will later set 
A = 1. This is justified because the short range single 
particle potential primarily serves to provide a boundary 
condition at the origin and its actual strength turns out 
to be unimportant. Next we expand the Greens function 
in a Taylor series in both 7 and A. Upon reorganizing the 
Taylor series expansion of Gn into an exponential series, 
we can factor out from the series expansion the unper- 
turbed iV-particle Greens function, the exact single par- 
ticle Greens function Gi(i) (defined through Eq. (fJJ with 
7 = 0) and exact two particle Greens function G2(i,j) 
(defined through Eq. (fj| with A = 0). Keeping just these 
terms, and neglecting those terms remaining in the se- 
ries expansion results in our 'first order' approximation 
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Gl{t)^.Gl{i,jy 



(5) 



where the superscript refers to the unperturbed so- 
lution. These turn out to be the dominant terms for 
all values of x, y, and in particular they reproduce the 
non-analytic behavior aforementioned. Note that the ap- 
proximation is 'first order' in the sense that the terms ne- 
glected in the series are all of order j 2 or higher; all terms 
of order 7 are kept, and some of the terms are summed up 
to infinite order in 7. The approximation turns out to be 
quite good even close to the Anderson transition where 
7 remains smaller than unity but clearly it can be im- 
proved systematically by repeating the above procedure 
for higher orders in 7. 

In order to obtain explicit results, we evaluate G\{i) 
and G 2 (i,j) in the limit 7 <C 1: 



Gi(x;t\y) 
Gi(x;t\y) 
G 2 (xi,x 2 ;t\yi,y 2 ) 
G 2 ) (x 1 ,x 2 ,t; \yi,y 2 ) 



w xr(2y) 

w T 7 (xi + x 2 )T y (\xi - x 2 \) 
x r( yi + y 2 )T(\ Vl - y 2 \) (6) 
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FIG. 1: (Color online) Density of the transmission eigenvalues 
for two values of disorder, V — 0.014 (solid red line) and 
r = 0.0285 (dashed blue line). The exponential gap at the 
origin increases with increasing disorder (decreasing F). 



where t(z) = tanh 2 z and 
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Here erfc'(x) = (y/Tr/2)e* erfc(x) where erfc(x) is the 
complementary error function. As a first check, we have 
verified that our solution for the distribution of the trans- 
mission levels obtained from Eq. ^ agrees with exact 
solutions known in the Q1D limit— for the special values 
of 7 = 1,2 and 4, for both small and large x where simple 
analytic expressions are available. In order to solve for 
the distribution of conductances using Eq. we con- 
sider the free-energy of N interacting particles on a line. 
First of all, we obtain the density of the eigenvalues <j(x) 
as a function of the disorder parameter T: 



da' 



<tV 2 7 ct ' ~ 4rin[cr'e/4] 



(8) 



The small and large x behavior are given by u{x) 
z£_ e -r'(x-x 3p ) f or x ^ x ^ anc j a ^ ~ 2LL fo r x ^> Xgp , 

Here x sp w [l/2r + 1] and T' w [T - 2r 2 ]. Note that 
r <C 1 in the insulating regime. Figure 1 shows the den- 
sity according to Eq. ([8|) for two different values of T, 
with exponential gap at the origin that increases with in- 
creasing disorder. In contrast, the density in the metallic 
regime is also constant, but starting at the originii. Thus 
our result suggests that the opening of a gap in the eigen- 
value spectrum could be considered as a signature of the 
metal-insulator Anderson transition. 



FIG. 2: (Color online) P(lng) in the insulating regime for 
three values of disorder, T = 0.014 (solid red line), T = 0.0285 
(dashed blue line) and F = 0.051 (dotted green line) corre- 
sponding to (lnc?) = -39.4, (lng) = -15.8 and (In g) = -8.9, 
respectively. The numerical data points are from Ref [2l]], for 
the same values of (lngr). 



The conductance in the strongly disordered regime 
7 1 is dominated by the smallest eigenvalue x\, al- 
though in contrast to Q1D it remains highly interacting. 
In this regime the distribution can be expressed in a sim- 
ple form: 



P(lng) = exp 

f(x) = v(x) 



-f 
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8e7Vr' 



=crfc 



(x-x sp )VF (9) 



where v(x) = [Fx 2 — In x — i In sinh 2:r] . Figure 2 shows 
plots for P(ln g) obtained from Eq. ^ for three values of 
disorder, compared with numerical data obtained from 
solving the tight binding Anderson model. In order to 
compare the results in more detail, we plot in Figure 3 
the variance a 2 of ln<?, which is consistent with a power 
law 2/3rd 7 . Figure 4 shows the third cumulant K3 which 
behaves as K3 ~ (— ln<?), and in the inset we plot the 
corresponding skewness \ = K3/ a^ 2 which seems to sat- 
urate at a value ~ 0.8. Note that the skewness changes 
sign around (— ln g) w 5, which is still far from the An- 
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derson transition which occurs^ around (— ln<; c ) w 1.3. 
Thus the distribution changes its shape from a large pos- 
itive asymmetry to a negative one that increases as the 
disorder is decreased toward the critical point. 



the properties of this matrix numerically allowed us to 
consider a simplified model with two independent ma- 
trix elements only. The distribution is then obtained 
as a function of a single disorder parameter that fixes 
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FIG. 3: variance of Ing plotted as a function of (— lngr) 2 ' 3 . 
The points are calculated from Eq. and the line is a fit. 

How close to the critical point can we approach with 
our present formulation? It turns out that once the skew- 
ness becomes negative, the contributions to the conduc- 
tance from eigenvalues other than the smallest one be- 
come important. Nevertheless, the value of the variance 
at the critical point (—lng) is 1.64, close to 1.09 ob- 
tained from numerical data£. This suggests that it should 
be possible to improve the calculations systematically to 
study the conductance distribution at the critical point 
by increasing the number of eigenvalues included in the 
calculation of the conductance distribution. 

In summary, we have developed a perturbative ap- 
proach that allowed us to obtain, from a transfer ma- 
trix formulation, the full distribution of conductances 
in the insulating regime of a 3D disordered conductor. 
The solution takes into account the interaction between 
the transmission eigenvalues that had been ignored in 
the past. The formulation involves a phenomenological 
matrix characterizing eigenvector correlations; analyzing 
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FIG. 4: The third cumulant K3 as a function of {— Ing). The 
inset shows the corresponding skewness \. Note that \ ( or 
K3) changes sign at {—In g) m 5. The Anderson transition 
occurs at (—Ing) ~ 1.3. 

the mean value (Ing), even though PQng) changes its 
shape from a positive to a negative skewness as the dis- 
order is decreased from the deep insulating regime toward 
the critical Anderson transition point. The results agree 
with recent numerical simulations of the full distribution. 
With current fabrication technology, it should be possi- 
ble to verify the predictions experimentally. While the 
method is developed for the insulating phase only, it also 
leads to a possible characterization of the Anderson tran- 
sition in terms of the opening of a gap in the spectrum 
of the transmission eigenvalues. By systematically im- 
proving the approximations, the method could therefore 
be used to study the qualitative features of the critical 
distribution near the Anderson transition point. 

We gratefully acknowledge helpful discussions with F. 
Evers and P. Markos. 
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